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Abstract — We analyze an uplink of a fast flat fading MIMO- 
CDMA channel in the case where the data symbol vector for each 
user follows an arbitrary distribution. The spectral efficiency of 
the channel with CSI at the receiver is evaluated analytically 
with the replica method. The main result is that the hierarchical 
decoupling principle holds in the MIMO-CDMA channel, i.e., 
the MIMO-CDMA channel is decoupled into a bank of single- 
user MIMO channels in the many-user limit, and each single- 
user MIMO channel is further decoupled into a bank of scalar 
Gaussian channels in the many-antenna limit for a fading model 
with a limited number of scatterers. 

I. Introduction 

Direct-sequence spread-spectrum code-division multiple ac- 
cess (CDMA) has been utilized as a multiple access scheme in 
wireless communication. As a method of overcoming the ca- 
pacity bottleneck in future wireless communication, multiple- 
input multiple-output (MIMO) systems have attracted so much 
attention since the latter half of the 1990s [1]. Recently, MIMO 
systems with CDMA technology (MIMO-CDMA) have been 
studied [2], [3]. 

Mantravadi et al. [2] evaluated the asymptotic spectral 
efficiency of a MIMO-CDMA channel with Gaussian modu- 
lation using the random matrix theory. However, other data 
modulation schemes such as quadrature phase shift keying 
(QPSK) modulation, or more generally, Af -quadrature am- 
plitude modulation (M-QAM) are commonly employed in 
practice. Thus, it is important to analyze the MIMO-CDMA 
channel with a modulation of the kind. The purpose of our 
study is to evaluate the spectral efficiency of the MIMO- 
CDMA channel in the case that the data symbol vector for 
each user follows an arbitrary distribution. 

It was reported [4] that the spectral efficiency of a CDMA 
channel with binary phase shift keying modulation can be 
evaluated with the so-called replica method. Then, the method 
was applied to the performance evaluation of MIMO systems 
[5], [6], [7]. Guo et al. [8] analyzed a CDMA channel in the 
case that the data symbol for each user follows an arbitrary 
distribution and claimed that the decoupling principle, the 
equivalent result to which has already been proved in the 
case of Gaussian modulation [9], [10], holds in an asymptotic 
limit, i.e., the CDMA channel is decoupled into a bank 
of scalar Gaussian channels. This principle makes possible 
the analytical evaluation of the spectral efficiency since the 



degree of freedom drastically decreases. However, it is not still 
clear except special cases [2] whether or not the decoupling 
principle holds in the MIMO-CDMA channel. 

In this paper, we claim that the decoupling principle holds in 
the MIMO-CDMA channel and evaluate the spectral efficiency 
of the channel using the replica method. 

II. Model 

We consider the uplink of a synchronous If -user MIMO- 
CDMA flat fading channel [2] 
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where the fcth user has M/. transmit antennas and the receiver 
has N receive antennas. x k — (x k , . . . ,x k VIk ) T is the data 
symbol vector for the fcth user and {s k ;l = 1,...,L} is 
the spreading sequence for the fcth user. We assume that 
{x k ; k — 1, . . . , K} are mutually independent complex ran- 
dom variables and that the real part and the imaginary part of 
{s k } are independent and identically distributed (i.i.d.) zero- 
mean random variables with the variance 1/2L. H k represents 
the N-by-Mh channel matrix for the fcth user, i.e., the (n, 
m/c)-element h k lmk of H k is the channel gain from the TOfcth 
transmit antenna of the fcth user to the nth receive antenna. 
We assume that {H k ; k — 1, . . . , K} are mutually indepen- 
dent. We consider the case that the noise is additive white 
Gaussian noise (AWGN), i.e., {nfj — 1,...,L} are i.i.d. 
zero-mean circularly symmetric complex Gaussian random 
variables with the covariance matrix a 2 I n, which are denoted 
by CJV(0, a 2 In)- Vi represents the iV-dimensional received 
signal vector. 

We write the entire data symbol vector x and the received 
signal vector y in a symbol period as x = (x lT , . . . , x K ) T , 



and y = (y^ 



y L ) , respectively. The maximal sum 



rate in the fast fading channel with perfect channel side 
information at the receiver is given by the conditional mutual 
information (the base of logarithms is taken to 2 in this paper) 
between x and y conditioned on the spreading sequences and 
the channel matrices [1] 
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where S = {sf ; I = 1, . . . , L, k = 1, . . . , if}, H = {H k ; k 
1, . . . , K}, and where p(y\x, S, Ti.) is given by 
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We define the spectral efficiency Cmimo-cdma as the maxi- 
mal sum rate per chip and per transmit antenna. In the many- 
user limit, where the number of users K and the spreading 
factor L tend to infinity with their ratio (3 = K/L fixed, the 
spectral efficiency is given by 
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where M = limfc^oa K ~ 1 X^fcLi Mf- and ft is the ratio of the 
average number of transmit antennas to the number of receive 
antennas, i.e., ft = M/N. The free energy T is defined as 



T=- lim —E[\rxp{y\S,U)\U\. 



(5) 



III. Replica Analysis of MIMO-CDMA 

We explain briefly the calculation procedure of the free 
energy (O [4], [8]. Substituting the following identity: 

d 
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to (0, we obtain 
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We assume that the limit with respect to K, L and the 
operation with respect to u are interchangeable. Then, (0 
becomes 
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Further, we assume that the result for positive integers u is 
valid for real number u in the calculation of ((9). 
We define A^-dimensional vectors v a as 
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where cc fe ' and # 
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are the original data 
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where X is defined as {x ,Q ; fc = 0, . . . , K , a = 



'} 



and S repserents {s k ; k = 1, . . . , K}. When K and L are 
sufficiently large with their ratio fixed, due to the central 
limit theorem, v = (v° , . . . 1 v uT ) T conditioned on X and 
H, follows approximately the zero-mean circularly symmetric 
complex Gaussian distribution with the covariance matrix 
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where w k is defined as [{H k x k -°) T , {H k x k < u ) T ] T . Then, 
we can evaluate (fTTT i and (fL2l as 

~P = ^lnE{e K ^ lGM ^}+0(K- 1 ), (14) 

G (tl) (Q) = -lndet(J + SQ) - A^uln^cr 2 ) - JVln(l + 
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where £ is defined as 
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where e„ is the w-dimensional vector whose elements are all 
one and (gi represents the Kronecker product. 

With a Hermitian matrix Q we define the moment generat- 
ing function of the data symbols of the fcth user as 

M ( k u) (Q) = E^* {exp [tr (QwV*)] } , (17) 

where X k represents {xfy*; mk = 1, . . . , Mk, a = 0, . . . , u}. 
Since Q satisfies the large deviation principle, with the saddle 
point method, ( [T4l is evaluated as 



lim = sup |~/3 -1 G (u) (Q) - 7 (u) (Q) 
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where the rate function is given by 
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Differentiating ( fT8l and ( fT9l with respect to Q and Q, 
respectively, we obtain the following equations giving extrema 
of CEH} and Q9): 



symbol vector and the replicated data symbol vector for the 
fcth user, respectively. {x k ' a , a — 0, ... ,u} are i.i.d. random 
variables following p(x k ) and {s k ;k — 1,...,K} are i.i.d. Q s 
random variables following p(s k ). Then, the expectation of 
[p(y\S,Ti.)] u with respect to y and S is given by 
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Differentiating (TT~8b with respect to u and substituting 
it, we can evaluate (|9]l as 



T = — lim 
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To evaluate the solution of d2"0l and (|2TT i in the u — > limit 
analytically we assume that the replica symmetry holds, i.e., 



Q s and Q are invariant under exchange of non-zero replica 
indexes. Then, Q s and Q can be written as 

Qs=( Q ° e -® M \ (23) 

\e u ®M* I u ®{Q l -Q) + e u el®Q) ' k ; 

Q S =( -A®* (24) 

\e u <E)M I U ®(Q -Q) + e u el®Q) 

where M, M are 7V-by-7V matrices and Q°, Q°, Q 1 , Q , 
Q, and Q are N-by-N Hermitian matrices. By solving < f20b 
we have in the u — ► limit 

Q° = 0, M = #\ Q 1 = Q - M, Q = R 1 R R 1 , 

(25) 

where Rq, R are defined as 

Ra = a 2 I N +f3(Q°-M-M* +Q), (26) 
R = a 2 I N + [3(Q l - Q). (27) 

It is straightforward to confirm that Rq, R are positive definite. 

We move on to calculating the moment generating func- 
tion A3. Ro can be decomposed into the product of two 
nonsingular matrices, i.e., Rq = ^/Rq-s/Rq . From ( |25] | the 
moment generating function (T% is evaluated as 
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where R a = R for a = 1, . . . , u 9 and where b k is defined as 

u 
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By applying the transform: 
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where p(y \x , H ;R) is defined as 
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Assuming that Rq is equal to R, from (l2Tb we obtain 
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Note that £ fe is averaged with respect to H k due to the law of 
large numbers. Substituting d33l to d26l i and ( f2Tb . we obtain 
the fixed-point equation 

i A ' 

/? = a 2 /jv+/3 lim — Y£ k , (36) 



k=l 



which is the extension of the Tse-Hanly equation [9] to the 
case of the MIMO-CDMA channel. 

Calculating G^{Q S ), I^iQ 5 ) and differentiating them 
with respect to u, we can evaluate T as 



In 2 



T 



lim 



K- 



1 K 

oo A ' — ' 

fe=l 



A r log(7rcr 2 e) 



+KL (CAT (0, a 2 /at) | |CAA (0, R)) , (37) 



where KL(-||-) represents the Kullback-Leibler divergence, 
and where C^ IMOl is defined as 
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Note that d37l i depends only on R. In the case that there 
exist multiple solutions of d36l l. one should choose the solution 
achieving the supremum of ( TT~8b in a neighborhood of u = 0, 
i.e., the solution minimizing (P37l i. 

From the above mentioned analysis we claim that the 
spectral efficiency of the MIMO-CDMA channel is evaluated 
as 

(-MIMO-CDMA — lim — > - l^fffc |C]\ 
K—>oo J\ — ' 
fc=l 
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EjffcpMiMoJ can ^ e interpreted as the spectral efficiency 
of the following single-user MIMO Gaussian channel for the 
fcth user: 

y k = H k x k + n k , n k ~ CW(0, R). (40) 

From the replica analysis for the moment sequence of the pos- 
5 ~r mean estimator (x^Jmimo-cdma = E[x£,Jj7, S, H], 
indicate that the MIMO-CDMA channel with the MMSE 
detector front end is decoupled into a bank of single-user 
MIMO channels with the MMSE detector front ends in the 
many-user limit. It is easy to extend this decoupling result to 
the cases of the linear MMSE detector or the matched filter. 

IV. Replica Analysis of MIMO 

So far we have not specified statistics of the elements of 
H k . In order to obtain a more concrete expression of d39l , we 
consider a fading model with a limited number of scatterers 
[6] 

H k = $ k A k ® k , (41) 

where ® k is a S k -by-Mk steering matrix which describes 
the propagation from the transmit antennas of the fcth user 



to Sk scattering objects between the fcth user and the receiver, 
where A k — diag(^, . . . , A% ) is a Sk-by-Sk diagonal 
matrix which accounts for attenuation at the scattering objects 
between the fcth user and the receiver, and where & k is a Sk- 
by-N steering matrix which describes the propagation from the 
scattering objects between the fcth user and the receiver to the 
receive antennas of the receiver. We assume that the elements 
of * A and ® k are i.i.d. random variables with the variances 
l/N and l/Sk, respectively, and that the elements of A k 
are independent random variables subject to the normalization 
E[tr(A k A k *)} = S k . 

It is difficult to evaluate (134-b and (1381 analytically except in 
special cases, e.g., the data symbol vector of the fcth user x k 
follows a circularly symmetric complex Gaussian distribution 
or the users and the receiver have a few numbers of the 
antennas. Thus, we evaluate ( |34l and ( |38l ) using the replica 
method in the many-antenna limit where Mk, Sk, and N tend 
to infinity with their ratios pk = Sk/N, jk = Mk/Sk fixed. 
One might think that the assumption of the many-antenna limit 
is impractical, but it can be a good approximate approach to 
systems with a few antennas if the elements of x k and H k 
follow circularly symmetric complex Gaussian distributions 
[1]. On the other hand, the spectral efficiency is not invariant 
under exchange of the order of the many-user limit and the 
many-antenna limit. In the MIMO-CDMA literature, however, 
it may be reasonable to take the many-user limit first. 

By regarding <l> fc and A k ® k x k as the channel matrix 
and the data symbol vector, respectively, we can evaluate the 
expectation of d38l > as 



limEjjfc [CmimoJ — _ 1™ E A fc@fc [Cmimo 2 ] 

b k ,M k ^oo 
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where pk — M^/N, where lim represents the many-antenna 
limit, and where Cmimo 2 * s me s P ectra l efficiency of the 
following MIMO channel for the fcth user: 



y k = A k @ k x k + h k , h k ~ CAT(0, C Is J. (43) 
£ fc2 and W k satisfy the following fixed-point equations: 
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mimo 2 = E[x k \y ,A ® } and where the op- 
erator E in (|45]> represents the expectation with respect to 
p(y k ,x k , A k ® k ). In the case that there exist multiple so- 
lutions of ( f44l > and (@5), one should choose the solution 
minimizing the spectral efficiency (1421 1. 

To evaluate (f34t we define a quantity -T^imo as 



f k 

-> MIMOi 
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where F is defined as 



F(X", H k ) = S k H k (x k - x k ^){x k 



x 



k:2 



Then, we obtain in the many-antenna limit [4] 
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Calculating the right-hand side of j49l , we can evaluate d34l > 

as 



km e k 
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where W is a solution of d44i > and d43b . We have added a 



~ k 

tilde to W K in order to make clear that W need not be equal 
to W k in (l42b . In the case that there exist multiple W one 
should choose the solution minimizing d42b . 

Furthermore, we can evaluate E[C^jimo ] 
applying the above mentioned method again. On the assump- 
tion that {x k lk ;rrik = 1, . . . , Mk} are mutually independent, 
E[Cmimo 2 1 ' n ^ C man y" antenna kmit is given by [6], [7] 
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where Cawgn is the spectral efficiency of the scalar Gaussian 
channel 
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£ fc2 and a k satisfy the following fixed-point equations: 
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where (a;*,.) awgn = E [x k - '' 



(55) 

mfc /AWGN = upmj |2/mJ- In the case that there exist 
multiple solutions of d53l . d541 l. and d55l l. one should choose 
the solution minimizing the spectral efficiency ( BTI ). 

The result of a replica analysis claims that the moment 
sequence of (xJ^^mimo-Cdma converges to the moment 
sequence of (i^Jawgn in the many-user and many-antenna 
limits. Since the MMSE detector is information lossless in 
the scalar Gaussian channel, Cawgn can interpreted as 
the spectral efficiency of the MMSE detector in the MIMO- 
CDMA channel. 

On the other hand, d45l > is evaluated as 



l s k 



ok 

c MIM0 2 „— o 



e +ik£% 



AWGNI^s fc l J 

(56) 



where Q WGN is a solution of (l53l l. d54l i. and (T55l l. In the case 
that there exist multiple f awgn one should choose the solution 
minimizing (l5TT i. We can obtain the same fixed-point equations 
for ( k and £ fc by means of calculating the expectation of 
Cmimo with respect to fc and then evaluating the asymptotic 
distribution of singular values of <& k . 

We can easily confirm that the spectral efficiency (|39l 
coincides with the spectral efficiency of the i.i.d. Rayleigh 
fading MIMO-CDMA channel in the case that scattering is 
very rich. In the 7^ — > limit with = pkjk fixed, (|53| > 
reduces to £ fc = ( k . Expanding d56l l with respect to 7fc and 
substituting it to d44l >. we obtain 

MS;£awgn(C & )In- 
2 (5?) 

(|57] > coincides with the fixed-point equation for ( k and W^ fe 
in the case of the i.i.d. Rayleigh fading. Hence, d39l converges 
to the spectral efficiency of the i.i.d. Rayleigh fading MIMO- 
CDMA channel in the 7 fc — > limit. 

The above mentioned results indicate that the MIMO chan- 
nel d40b with the MMSE detector front end is decoupled into 
the bank of scalar Gaussian channels with the MMSE detector 
front ends in the many-antenna limit even if the number of the 
scatterers is limited. It is easy to extend this decoupling result 
to the cases of the linear MMSE detector or the matched filter. 

V. Numerical evaluation 



( k 2 = lim 



N 



-tr 



W = R 



To evaluate the spectral efficiency numerically we consider 
the case that p = pk, 7 = 7&, stre i.i.d. random 

variables with the variance P = E[|a;^J 2 ], and 1^4^ | 2 = 1 
with probability 1. p = p-f represents the ratio of the number 
of transmit antennas to the number of receive antennas, i.e., 
/i = M/N. We denote the received signal-to-noise ratio per 
transmit antenna by SNR = P/(Jq. Figure Q] displays the 
spectral efficiency d39i > versus SNR for QPSK modulation (4- 
QAM). We find that the spectral efficiency of the MIMO- 
CDMA channel with QPSK modulation is very close to the 
capacity but there exists a large gap between the spectral 
efficiency of the MMSE detector and the capacity in the highly 
loaded system. The degradation of the spectral efficiency due 
to the decrease of scatterers is shown in Fig. [2] In this case, an 
interesting observation is that the spectral efficiency does not 
degrade so much compared with the rich scattered environment 
(7 — ► 0) even when the number of scatterers is comparable 
with the number of the antennas. 

VI. Conclusion 

We evaluated the spectral efficiency of the uplink of the 
MIMO-CDMA channel using the replica method. The main 
result is that the hierarchical decoupling principle holds in 
the MIMO-CDMA channel, i.e., the MIMO-CDMA channel 
is decoupled into the bank of single-user MIMO channels in 
the many-user limit. The resulting single-user MIMO channel 
is further decoupled into the bank of scalar Gaussian channels 
in the many-antenna limit for the fading model with a limited 
number of scatterers. We found numerically that the spectral 



efficiency of the MIMO-CDMA channel with QPSK modula- 
tion is very close to the capacity but there exists a large gap 
between the spectral efficiency of the MMSE detector and the 
capacity in a highly loaded system. 
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Fig. 1. Spectral efficiency of the MIMO-CDMA channel versus SNR. 
ii = 1.0 and 7 = 1.0. 
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Fig. 2. Spectral efficiency of the MIMO-CDMA channel for QPSK 
modulation versus 7. fj, = 1.0 and SNR = 10 dB. The spectral efficiency 
coincides with the spectral efficiency of the i.i.d. Rayleigh fading MIMO- 
CDMA channel in the 7 — * limit. 



